Quantum Billiards with Surface Scattering: Ballistic Sigma-Model Approach 
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Statistical properties of energy levels and eigenfunctions in a ballistic system with diffusive 
surface scattering are investigated. The two-level correlation function, the level number variance, 
the correlation function of wavefunction intensities, and the inverse participation ratio are calculated. 
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The statistical properties of spectra of disordered dif- 
fusive systems are now well understood. Using the su- 
persymmetric cr-model approach it has been possible to 
demonstrate the relevance of the random matrix theory 
(RMT) and to calculate deviations from its predictions 
both for the level and eigenfunction ^-|^ statistics. 
Generalization of these results to the case of a chaotic bal- 
listic system (i.e. quantum billiard) has become a topic of 
great research interest. For ballistic disordered systems 
the fj-model has been proposed 0| , with the Liouville op- 
erator replacing the diffusion operator in the action. It 
has also been conjectured that the same cr-model in the 
limit of vanishing disorder adequately describes statisti- 
cal properties of spectra of individual classically chaotic 
system. This conjecture was further developed in 
where the cr-model was obtained by means of energy av- 
eraging, the Liouville operator replaced by its regulariza- 
tion — the Perron-Frobenius operator and some neces- 
sary conditions for validity of this description were put 
forward. 

However, straightforward application of the results of 
Refs. lll-^j^ to the case of an individual chaotic system 
is complicated by the fact that the eigenvalues of the 
Perron-Frobenius operator are unknown, while its eigen- 
functions are extremely singular. For this reason the cr- 
model approach has so far failed to provide explicit re- 
sults for any particular ballistic system. 

To overcome this difficulty, we consider a model of a 
billiard with surface disorder leading to diffusive scat- 
tering of a particle in each collision with the boundary. 
This models behavior of a quantum particle in a box 
with a rough boundary which is irregular on the scale 
of the wave length. Since the particle loses memory of 
its direction of motion after a single collision, this model 
describes a limit of an "extremely chaotic" ballistic sys- 
tem, with the typical relaxation time being of order of the 
flight time. (This should be contrasted with the case of a 
relatively slight distortion of an integrable billiard JTof.) 
One might naively think that all results for such a model 
could be obtained by setting Z « L in a system with bulk 
disorder. In fact, the level statistics in a system with bulk 
disorder and arbitrary relation between mean free path I 



and system size L were studied in |ll| . We will see, how- 
ever, that our results are qualitatively different in some 
respects, which shows that systems with bulk and sur- 
face disorder are not equivalent. On the other hand, our 
findings are in agreement with general expectations for 
chaotic billiards based on a semiclassical (trace formula) 
treatment ]T^ . 

To simplify the calculations, we consider a circular bil- 
liard. A similar problem was studied numerically in Ref. 
JlSf for a square geometry. We consider only the case of 
unitary symmetry (broken time-reversal invariance) ; gen- 
eralization to the orthogonal case is straightforward and 
will be given elsewhere The level statistics for the 
same problem were independently studied in Ref. p5| . 

Properties of the Liouville operator. Our starting 
point is the sigma-model derived recently ^ for ballistic 
disordered systems. The effective action for this model 
has the form 
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Here g(r,n) is a 8 x 8 supermatrix, which depends on 
the coordinate r and direction of the momentum n. 
The angular braces denote averaging over n: {0{n)) = 
J dnO{n) with the normalization J dn — 1. The ma- 
trix g is constrained by the condition g{r,n)^ — 1, 
and can be represented as <? = UAU~^, with A = 
diag(l, 1, 1, 1, —1, —1, —1, —1); see [0^] for more detailed 
definitions. Since we are interested in the clean limit with 
no disorder in the bulk, the second term in the action (|l|) 
containing the elastic mean free time r is zero everywhere 
except at the boundary where it modifies the boundary 
condition (see below). 

Most of the statistical properties of energy levels 
Jl6|,p|J^ and eigenfunctions |^,^| are determined by the 
structure of the action in the vicinity of the homoge- 
neous configuration of the g-field, g{r,n) = A. Writing 
U = I — W/2 + . . ., one finds the action in the leading 
order in W to be 
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where the indices 1 , 2 refer to the "advanced-retarded" 
decomposition of W, and the Liouville operator K = 
vpnV . This "Unearized" action has the same form as 
that of a diffusive system, with the diffusion operator be- 
ing replaced by the LiouviUe operator. This enables us 
to use the results derived for the diffusive case by substi- 
tuting the eigenvalues and eigenfunctions of the operator 
K for those of the diffusion operator. 

The operator K should be supplemented by a bound- 
ary condition, which depends on the form of the surface 
roughness. As a model approximation we consider purely 
diffusive scattering for which the distribution func- 
tion iy9(r, n) of the outgoing particles is constant and is 
fixed by flux conservation: 

tp{r,n) =TT {Nn')ip{r,n')dn', (TVn) < 0. 

J{Nn')>0 

Here the point r lies at the surface, and N is an outward 
normal to the surface. This boundary condition should 
be satisfied by the eigenfunctions of K. 

The eigenvalues A of the operator K corresponding to 
angular momentum I obey the equation 
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d0 sin exp [2iW + 2^ sin 9] = 0, (3) 



where ^ = RX/vp, and R is the radius of the circle. For 
each value of / = 0, ±1,±2, . . . Eq.(||) has a set of so- 
lutions ^ik with ^ik = and ^a- = C-fc- These fa- 
can be naturally labeled with fc = 0, ±1, ±2, . . . (even I) 
or fc = ±1/2, ±3/2, . . . (odd /). For / = fc = we have 
Coo = 0, corresponding to the zero mode (p{r, n) = const. 
All other eigenvalues have positive real part Re fa > 
and govern the relaxation of the corresponding classi- 
cal system to the homogeneous distribution in the phase 
space. 
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FIG. 1. First 11 x 11 (0 < fc,/ < 11) eigenvalues of the 
Liouville operator K in units of the inverse time of flight, 
vp/R. 

The asymptotic form of the solutions of Eq. (||) for large 



|fc| and/or |Z| can be obtained by using the saddle-point 
method, 
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0.66; + 0.141n/ + 0.557rifc, 0<fc<Z 



(lnfc)/4 + 7ri(fc + 1/8), 



< / < fc 



(4) 



Note that for fc = all eigenvalues are real, while for 
high values of fc they lie close to the imaginary axis and 
do not depend on I (see Fig. 1). 

Level statistics, low frequencies. We define the 
level correlation function in a standard way, 

R2{uj) = {AVf{j^{e + ujHe))~l, 

where i'{e) is the density of states, A — {Vi')^^ is the 
mean level spacing and V = ttR^ is the area. In the 
range of relatively low frequencies (which for our prob- 
lem means uj ^ vp/R, see below) the function R2{oj) 
quite generally has the form B 



i?2(s) = 6{s) - {tts)- 



sm TTs 



A (i? A /nvp ) sin^ tts , 



(5) 



where s — uj/A. The first two terms correspond to the 
zero-mode approximation and are given by RMT, while 
the last one represents the non-universal correction to the 
RMT results. It contains information about the opera- 
tor K through the dimensionless constant A = X^'CaT;^' 
where the prime indicates that the eigenvalue foo = is 
excluded. The value of A, as well as the high-frequency 
behavior of the level correlator (see below), can be found 
from the Altshuler-Shklovskii type spectral function [Hm: 



Y^{\u-^oJ)-\ (6) 
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According to the Cauchy theorem, Si can be represented 
as an integral in the complex plane. 



'^'^^ \vp) 2m fc (z-iujR/vp)^ Ji{z) 



where the contour C encloses all zeroes of the function 
Ji{z). Evaluating the residue a,t z — iuR/vp, we find 



Siico) = -{R/vpf — 



\nJi{z). 



(7) 



z—iujR/ vf 



Considering the limit — > and subtracting the contri- 
bution of Aoo = 0, we get 

A = -19/27 - 1757rVll52 + 64/(97r2) « -1.48. (8) 

In contrast to the diffusive case, this constant is negative: 
the level repulsion is enhanced with respect to result for 
RMT. Eq.(|^) is valid as long as the correction is small 
compared to the RMT result, i.e. provided w is below 
the inverse time of flight, vp/R, which plays the role of 
the Thouless energy for our problem. 
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Level statistics, high frequencies. In the range 
uj ^ A the level correlation function can be decomposed 
into the smooth Altshuler-Shklovskii part = 
(AV27r2)ReS'(w) and the part i?f ° which oscillates 
on the scale of the level spacing. Evaluating the asymp- 
totic behavior of Si{uj) from Eq.(|7|), we find in the high- 
frequency regime when uj 3> vp/R- 



AR 

Vp 



,1/2 f UjR TT 
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\2ttujRJ \ Vp 4 



(9) 



The oscillating part of the level correlation function 
R2^'^{s) for frequencies uj ^ A is given by M 



R2^^{s) = {1/271^) cos{2tts)D{s), 



where D(s) is the spectral determinant, 



(10) 



Dis) = Yl {l-isA/X 



kl 



r\l + isA/\ki)-\ 



Since A~^9^ lnL'(s)/9s^ — — 2ReS'(w), we can restore 
D{s) from Eqs.(||), (|^), up to a factor of the form exp(ci-|- 
C2s), with ci and C2 being arbitrary constants. These 
constants are fixed by the requirement that Eq.(|o|) in 
the range A <^ lo <^ vp / R should reproduce the low- 
frequency behavior (^. As a result, we obtain 



^ ' V 2/ TV J- J- J, 



Ji(isiV-i/2)j;(_isAr-i/2) 



(11) 



TT^SzCs) = 1 + 7 + ln 
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when s > iV^/^. Here iV = [vp/RAf = {ppR/2)^ is the 
number of electrons below the Fermi level, 7 ~ 0.577 is 
Euler's constant, and A is defined by Eq.(||). The first 
three terms at the rhs of Eq. (|lj) represent the RMT con- 
tribution (curve 1 in Fig. 2). The two asymptotics (|lj) 
and ( p^ are shown in Fig. 1 as curves 2 and 3 respec- 
tively. 
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FIG. 2. Level number variance T.2{E) as a function of en- 
ergy; s = E/ A. Curve 1 shows the RMT result, while curves 
2 and 3 correspond to asymptotic regimes of low ( p^ ) and 
high (^^ frequencies. The saturation value Ej"' is given in 
the text. 



For high frequencies to ^ vp/R this yields the following 
expression for the oscillating part of the level correlation 
function: 



TT* /ARV /2nuj 
= ens 



i?0^'=(tj) = — I I COS , ^ 

^ ^ ^ 128 \ VP J V A 



(12) 



It is remarkable that the amplitude of the oscillating 
part does not depend on frequency. This is in contrast 
to the diffusive case, where in the Altshuler-Shklovskii 
regime {ui above the Thouless energy) the oscillating part 
R2^'^{uj) is exponentially small ||]. 

The level number variance. The smooth part of 
the level correlation function can be best illustrated by 
plotting the variance of the number of levels in an energy 
interval of width E = sA, 



(13) 



S2(s) = / (s- |S|)i?2(S)dS, 



A direct calculation gives the following asymptotic be- 
havior: 

t:^^2{s) = 1 + 7 + ln(27rs) + As^/{2N) (14) 
when s A^-^/^ and 



As seen from Fig. 2, the two curves (14) and ( |15| ) per- 
fectly match in the intermediate regime, s ~ N^, and 
taken together they provide a complete description of 
E2(s). According to Eq.(p^, the level number variance 
saturates at the value S^"^ = TT-^{l+j+ln{16N^^'^ /n'^J 
in contrast to the behavior found for diffusive systems | 
or ballistic systems with weak bulk disorder ||ll[ . The sat- 
uration occurs at energies s A^^/^, or in conventional 
units E ^ vp/R. This behavior of £2(3) is expected for 
a generic chaotic billiard |lj]. It is also in good agree- 
ment with the results for ^2(5) found numerically for a 
tight-binding model with moderately strong disorder on 
boundary sites ]l3t . 

Eigenfunction statistics. Now we study correlations 
of the amplitudes of an eigenfunction in two different 
points. Correlation of different eigcnfunctions will be 
considered elsewhere [Q . Following Ref . ^ , we define 

a(ri,r2,i?) = AV^i^, \Mri)^M\^SiE - e^)), 

where i/'/i ^tre the eigcnfunctions corresponding to the ex- 
act single-particle states /i. A calculation analogous to 
that of Ref. [D yields 



a{ri,r2,E) = 1 4- n(ri, ra). 



(16) 
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where 11 is the Green's function of the operator K 
integrated over directions of momentum, n(T"i,r2) — 
J dnidn2 g{ri,ni;r2,n2). Here g is the full Green's 
function of the operator K, i.e. a solution to the equa- 
tion 



(17) 



Kg{ri,nx;r2,n2) = 
{t::^)-' [S{ri - r2)(5(ni - na) - V~^] . 

Direct calculation gives: 

n(ri,r2) =ni(ri,r2)+n2(ri,r2), 

Tli{ri,r2) =k4ri-r2)-V-^ [ dr[~k4r[ ~ r2) (18) 



n2(ri,r2) 
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E 

k=l 



4fc2 



1 / ?'i rn \ ^ 



V i?2 J 



COS fc ( 



where kd{r) = l/{TTpp\r\), and (r, 6*) are the polar co- 
ordinates. This formula has a clear interpretation. The 
function 11 can be represented as a sum over all possi- 
ble paths leading from ri to r2, with possible surface 
scattering in between. In particular, the function Hi 
corresponds to trajectories coming directly from ri to 
r2 with no reflection from the surface. Therefore, the 
term Hi is universal and is not sensitive to the geometry 
of the system. It can be obtained from the RMT-like 
prediction that amplitudes of different wavefunctions are 
independent Gaussian variables |jl^. More precisely, we 
find that the function kci{ri — r2) = Jq{pf\tci — i"2|) of 
Refs. j^Jl^ is replaced in Eq.(p^) by its smoothed ver- 
sion, kd{ri — r2) — 1/{ttpf\j^i — r2|)- This is because 
our semiclassical approach is valid on scales much larger 
than the wave length. To cure this flaw, one has to re- 
place fcd(r) by kd{r) in the expression Eq.(|lJ) for Hi. 
The second term, 112, is due to the surface scattering. 
It can be shown that in the numerator 4fc2 — 1 the 
first term comes from trajectories with only one surface 
reflection, while the second sums up contributions from 
multiple reflections. A formula analogous to (^6|) was 
proposed very recently for a generic chaotic system ||l9[| . 

Finally, we calculate the inverse participation ratio, 
(P2) = V^^ j dra{r,r), which characterizes the degree 
of spatial uniformity of eigenfunctions. The RMT pre- 
diction for this quantity, Pj^"* = 2/V is recovered from 
Eqs.(^^, (|8|) if we take into account the first term in 
the expression for Hi, since kd{0) = 1. The leading cor- 
rection comes from the single-reflection contribution to 
the term 112, and is equal to 



SP2 = V-\4TrpFRy^HpFR) 



In conclusion, we have used the ballistic cr-model ap- 
proach to study statistical properties of levels and eigen- 
functions in a billiard with diffusive surface scattering, 



which exemplifles a ballistic system in the regime of 
strong chaos. We have found that the level repulsion and 
the spectral rigidity are enhanced compared to RMT. 
In particular, the level number variance saturates at the 
scale of the inverse time of flight, in agreement with 
Berry's prediction for a generic chaotic system jl^. As 
another manifestation of the strong spectral rigidity, the 
oscillating part of the level correlation function does not 
vanish at large level separation. We have also consid- 
ered correlations of eigenfunction amplitudes in different 
spatial points and calculated the deviation of the inverse 
participation ratio from the RMT value. 
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